Conical Emission from Shock Waves in Ne(l— 20 AGeV)+U Collisions 
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The formation and propagation of high-density compression waves, e.g. Mach shock waves, in 
cold nuclear matter is studied by simulating high-energy nucleus-nucleus collisions of Ne with U 
in the energy range from i5iab = 0.5 AGeV to 20 AGeV. In an ideal hydrodynamic approach, the 
high-density shock wave created by the small Ne nucleus passing through the heavy U nucleus is 
followed by a slower and more dilute Mach shock wave which causes conical emission of particles 
at the Mach cone angle. The conical emission originates from low-density regions with a small flow 
velocity comparable to the speed of sound. Moreover, it is shown that the angular distributions of 
emitted baryons clearly distinguish between a hydrodynamic approach and binary cascade processes 
used in the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport model. 



in 



> 
m 

rn 
O 

o 



X 



I. INTRODUCTION 

Hydrodynamic models predict a sideward emission of 
nuclear matter in fast nucleus-nucleus collisions due to 
the transformation of kinetic energy of the projectile 
into compression and heat energy of the medium. High- 
density shock waves created during the collision of un- 
equal nuclei push matter in transverse direction, gener- 
ating a measurable preferential emission at a well-defined 
Mach angle [iHl] ■ This angle is connected to the medium 
properties, in particular to the velocity of sound 



cl = dp/de, 
by the classical Mach cone formula 



(1) 



(2) 



Here, Ush denotes the velocity of the leading head shock 
wave traveling through the target nucleus. 

In recent years, the investigation of (Mach) shock 
waves has re gain ed tremendous attention, both on the 
theoretical [7Hl9l[ as well as on the experimental side [20l - 
[28j . considering shock waves from partonic projectiles. 

In this paper, we return to the original idea and re- 
investigate the creation of Mach-likc shock waves as 
well as the sideward deflection of matter in asymmetric 
nucleus-nucleus collisions using (3-|-l)-dimensional ideal 
hydrodynamics [29j and the Ultra-relativistic Quantum 
Molecular Dynamics (UrQMD) transport model [13, ■ 
The conical emission at Mach angles 6'mc (where 9 de- 
notes the polar angle between the beam axis and the flux 
of matter) persists over a wide range of beam energies and 
impact parameters, allowing for a comprehensive study 
even at lower beam energies. We analyze the change of 



the particle emission angle 0iab with increasing impact 
parameter and show the dependence of the Mach shock 
wave on the size of the projectile nucleus. 

The (3-1-1 )-dimensional ideal hydrodynamic calcula- 
tions are performed for central and non-central Ne-|-U 
collisions at 0.5-20 AGeV, applying a chiral hadronic 
equation of state (EoS) that exhibits a phase transition 
to a chirally restored phase, as it is expected by lattice 
calculations [32l - [3^ . Therefore, we are able to investi- 
gate the impact of this phase transition on the emission 
angle of particles from Mach-like shock waves in heavy 
ion collisions. 

Additionally, analogous calculations are performed us- 
ing the UrQMD transport model. The comparison of 
the different models with experimental data will allow 
to draw conclusions about the underlying process of the 
conical emission. While the conical emission of the reac- 
tion products in hydrodynamics is due to the formation of 
shock waves in the nuclear matter, the emission pattern 
in the UrQMD model is generated by binary collisions 
and does not show a Mach cone pattern. 



II. HYDRODYNAMICS 

Each hydrodynamic simulation of a collision process is 
initialized when the Ne nucleus starts to penetrate the U 
nucleus and is performed in the rest frame of the Uranium 
target. The ground state energy density distribution of 
each nucleus is given by the Woods-Saxon form 



= eo/ 



1 + e('-''°)/'* 



(3) 
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where eo = 147 MeV/fm'^ denotes the ground-state en- 
ergy density. 

The mean radius of a nucleus of mass A is computed 
via ro = (l.l2^(i/3) _ 0.86 ^("^Z^)) fj^ ^nd the surface 
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Figure 1: (Color online) The sound velocity Ca of the EoS 
used in our calculations as a function of e/eo and n/no. 



thickness is set to d — 0.54 fm [35|- In the longitudi- 
nal direction the density distribution of the projectile is 
Lorentz contracted. The distribution of the net baryon 
number density in the cold nuclei is computed accord- 
ingly with TiQ = 0.16 fm~^. 

Ideal hydrodynamics denotes the conservation of en- 
ergy, momentum, and net baryon number 



0, 



inu'^) = 0. 



(4) 



These equations are numerically solved by the SHASTA 
algorithm [2^ [36l-l39^ on a discretized grid with the ideal 
hydrodynamic energy-momentum tensor given by T^'' 



{e+p)i 



7(l,v),7 = (l-z;2) 



-pg'^'"'^ (with the four- velocity of the fluid 



'1/2 



and the net baryon density n). 



The algorithm ensures the flux-corrected transport of the 
thermodynamic quantities over the computational grid as 
described in Refs. (36l - [38| . The symmetric Eulerian grid 
is flxed in the target (computational) frame with a static 
cell size Ax = 0.2 fm. Time steps are set to Ai — 0.08 fm. 
Changing these numbers does not significantly affect the 
results. 

The pressure p is connected to the energy density e and 
the net baryon number density n by the EoS p(e, n). For 
this study a hadronic EoS, derived from a chiral hadronic 
SU(3) Lagrangian which includes the lowest baryon octet 
together with the multiplcts of scalar jjseudo-scalar, vec- 
tor, and axial- vector mesons, is used [40|. All parameters 
of the model EoS are fixed either by symmetry relations, 
hadronic vacuum observables or nuclear matter satura- 
tion properties. The model exhibits a smooth decrease 
of the chiral condensates (cross o yer) for high temper- 
ature and low baryonic potential [U In addition, 
the model also provides a satisfactory description of re- 
alistic (finite-size and isospin asymmetric) nuclei and of 
neutron stars [4l|, IH, Q ■ Additional baryonic degrees 
of freedom change the smooth cross over into a first- 
order phase transition with a critical end point (CEP) 
at T ~ 180 MeV, ~ 115 MeV, depending on the 
couplings dO, in, l45n47| . These values for the CEP 



are within the range expected from lattice QCD calcula- 
tions Is^-ls^- lisj . For a detailed discussion of this EoS see 
Refs. do, llj. The velocity of sound Cs in the n-e-plane 
of this EoS is depicted in Fig. [T] 

As the projectile nucleus hits the target with super- 
sonic velocity, the nuclear matter in the collision zone 
gets highly compressed and thereby creates a strong 
shock wave that moves with Vs\i > Cs through the tar- 
get nucleus. The shock velocity 



{P2 -pi){ei +P2) 



-,1/2 



(62 - ei)(ei +pi) 



(5) 



can be derived from the one-dimensional relativistic 
shock adiabat {Tauh adiabat [1^). Quantities with in- 
dex 1 denote the unperturbed nuclear matter while the 
index 2 denotes the energy and pressure in the compres- 
sion zone. Unlike the Glassgold-Heckroth- Watson ap- 
proach [531 and the participant-spectator model [51j, in 
head-on collisions the evolving violent shock wave com- 
pletely destroys the entire target nucleus, leaving no tar- 
get fragments in the target frame at all. The very hot, 
dense, and fast moving region near the collision axis is 
commonly referred to as head shock [l], |3, [^, [l6j . The 
center of the head shock reaches the following maximal 
local rest frame (Irf) energy and baryon number densities 
during the early evolution of the shock wave 

beam energy: Eu^ [AGeV] 0.5 1 10 20 



Cmax/eo (Irf) 

«max/«0 (Irf) 



3.5 13 32 50 
3 10 11 11. 



The further evolution is accompanied by a flow of energy 
to the outer regions of the system. This causes a signif- 
icant weakening of the head shock. For beam energies 
Eiah ^ 1 AGeV the chirally restored phase is reached 
in the head shock. At the back of the shock wave, a 
slowly propagating region (v ^ c^) with low energy den- 
sity is formed. In the head shock, matter is pushed ahead 
while it gets deflected sidewards in the outer regions of 
the shock wave, leading to a Mach shock structure in the 
medium 0, Q. 

Figure [2] shows the distribution of the laboratory frame 
energy density eiab in the reaction plane at different 
stages of a non-central (6 = 7 fm) collision of Ne-|-U at 
E\ah = 10 AGeV. The upper panel displays the collision 
process in hydrodynamics, the lower panel in the UrQMD 
framework. In hydrodynamics, the evolving shock wave 
is clearly visible. Note that the laboratory energy density 
is cut off at 7 GeV/fm'^ in this figure. 

The described shock wave creation is comparable to 
the scenario of a fast parton jet which distributes energy 
and momentum to the medium [7l-[l9l|. Particles emit- 
ted directly from the head shock have large momenta in 
forward direction and therefore are emitted at small 0iab 
angles, while particles emitted from outer regions of the 
(Mach) shock wave have low momenta and are predomi- 
nantly emitted in sideward directions. 
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Figure 2; (Color online) Laboratory energy density (eiab) distribution in the reaction plane at different stages of a Ne-|-U 
collision at _Eiab ~ 10 AGeV and b — 7 fm, calculated using the hydrodynamic (upper panel) and the UrQMD framework (lower 
panel) from initial to final state. Note that in this plot the laboratory energy density is cut off at 7 GeV/fm'^. 



The hydrodynamic evohition is stopped at time tf. 
Here, tf is chosen large enough to ensure the full transi- 
tion of the head shock through the target nucleus. The 
energy density distribution in the reaction plane in this 
final state </ is shown in Fig.[3]for central Ne+U collisions 
at different beam energies, calculated using the hydrody- 
namic approach. The shock wave becomes stronger with 
increasing beam energy and thus increasing energy den- 
sity in the head shock. In fact, the energy density reaches 
values up to eiab = 11 GeV/fm'^ at tf for the highest 
beam energy of Euh = 20 AGeV (cf. Fig. [S]). Thus, the 
velocity of the head shock comes close to Vsh = 1 accord- 
ing to Eq. (O. 

The decoupling of the hydrodynamic system into par- 
ticles is done at t/ performing an isochronous freeze-out. 
In a first scenario, we consider a free streaming freeze- 
out (FS FO). This is a rough approximation for the dis- 
tribution of light (He) and intermediate mass (Li, Be) 
reaction products which exhibit much less 1 /A) ther- 
mal smearing than nucleons 0, ij. The spectra of the 
reaction products are obtained by computing the rela- 
tivistic kinetic energy of the nucleons from every cell on 
the computational grid {i,j,k) via 



N 



(6) 



with a nucleon mass of = 939 MeV and the Lorentz 
gamma factor 



-1/2 



(7) 



where vl'^''^ denotes the velocity and n*'^'*^ the net baryon 
density of the specific fiuid element. The polar angle 



between the matter fiux and the beam axis 



COS 



(8) 



is computed for each cell and all baryons on the grid. I.e., 
thermal smearing and resonance decays are neglected in 
this scenario. 

A second freeze-out scenario is the Cooper-Frye pre- 
scription ^] (CF FO) 



dN 

E-^ = I da^p''f{x,p), 



(9) 



where f{x,p) denotes the boosted distribution function 
and da^ = {d^x, 0) the normal vector on the hypersur- 
face. We will consider the Cooper-Frye freeze-out with 
and without final state interactions and resonance de- 
cays which are performed using the UrQMD model as an 
afterburner. The details about performing this Cooper- 
Frye freeze-out after the hydrodynamic evolution as well 
as the integration into the UrQMD model are described 
in detail in Ref. (ssj . The hydrodynamic calculations 
are contrasted with simulations within a binary collision 
hadronic cascade, the default UrQMD model without 
a hydrodynamic phase [30l [3l| . In the latter two ap- 
proaches the emission angle of the baryons with respect 
to the beam axis is computed via ^lab = cos~^ {pz/\p\). 

The hydrodynamic calculations show a sizeable frac- 
tion of particles with low kinetic energies. They are emit- 
ted from regions with low flow velocities behind the lead- 
ing head shock at angles in the range of 9mc (Eq. (H])). 
In contrast, particles with high kinetic energies stemming 
from the highly compressed head shock should be emitted 
under small angles. Since the head shock gets remark- 
ably stronger with increasing beam energy, an increasing 
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Figure 3: (Color online) Laboratory frame energy density 
eiab [GeV/fm^] distribution in the reaction plane at tf, i.e. 
when the head shock wave has passed through the tar- 
get nucleus, calculated in a hydrodynamic approach. The 
maximum energy density in the center of the head shock 
is eiab ~ 0.2, 1.6, 2, 11 GeV/fm^ for projectile energies of 
-Eiab = 0.5, 1, 10, 20 AGeV. 



number of high energetic particles is expected at small 
angles and higher beam energies. 



III. RESULTS AND DISCUSSION 

Figure ID shows the kinetic energy distributions of mi- 
cleons, obtained with the free streaming freeze-out, for 
different beam energies. Beam energies in the range of 
-Eiab < 10 AGeV show a distinct peak for nucleons with 
low kinetic energies (i?kin < 50 MeV). For higher beam 
energies, however, the peak height decreases, but due to 
the creation of a stronger head shock more particles are 
produced at higher i?kin- 

This strong head shock, leading to a particle emission 
in forward direction, can clearly be seen from the angu- 
lar distribution of the emitted nucleons in the uppermost 
panel of Fig. [SJ where the hydrodynamic scenario was 
applied, followed by the free streaming freeze-out. While 
one clearly observes the suppression of particle emission 
under small angles for £^iab < 10 AGeV, a strong head 
shock develops at higher beam energies and pushes a con- 
siderable amount of matter ahead, resulting in a particle 
emission at small angles {Ouh < 20°) and a reduced emis- 
sion of particles in ^MC-direction. However, for all beam 
energies investigated in this study, a considerable number 
of nucleons is emitted at large angles (50° < 6'iab < 80°) 
due to the sideward deflection of matter in the evolving 
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Figure 4: (Color online) Energy spectra of reaction prod- 
ucts in Ne-|-U collisions. The majority of nucleons at the 
Mach angle arises from the low-energy (shaded) region -Ekin < 
100 MeV, corresponding to a nucleon velocity v — 0.43. High- 
energy particles, however result from head shock regions. 



Mach shock wave. 

The two middle panels in Fig. [S] show results from hy- 
drodynamic calculations with particle production via the 
Cooper-Frye freeze-out (second panel from top) and sub- 
sequently the UrQMD afterburner (third panel). The 
maxima of the distributions ^max are now located at 
much smaller angles within a range 20° < 6'iab < 
50°. Here, the fully isotropic, but Lorentz-boostcd ther- 
mal particle distribution considered in the Cooper-Frye 
freeze-out (ssj is superimposed on the results from the 
hydrodynamic calculations. 

To understand this result one can consider a simple 
toy model where a significantly strong isotropic particle 
distribution, boosted in beam direction, is superimposed 
on the results from pure hydrodynamics. This simple 
approach can reproduce the effect of the Cooper-Frye 
freeze-out on the particle distribution. Note, that due 
to reasons of geometry, a fully isotropic (thermal) distri- 
bution has a sm(0)-shape in the conioving frame when 
plotted as dN/dO. If this distribution is moving with re- 
spect to the observer (laboratory frame) the sinusodial 
shape is deformed and the maximum of the distribution 
is shifted to smaller angles. 

Particle rescattering and the decay of baryonic reso- 
nances within the UrQMD afterburner, however, shift 
the maxima of the angular distribution of emitted nucle- 
ons to slightly larger angles (cf. third panel of Fig. [S])- 

The lowest panel of Fig. [5] shows the results of default 
UrQMD calculations. The momentum conserving parti- 
cle scattering results in much broader and smoother an- 
gular distributions, with maxima about 30° < ^lab < 70°, 
dependent on the beam energy. The shape of the curves 
simply reflects an isotropic distribution that is Lorentz- 
boosted with a certain velocity on top of the fermi dis- 
tributed nucleons which did not interact at all. 

Figure [5] illustrates that a cut in the kinetic energy 
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Figure 5: (Color online) Angular distributions of nucleons 
obtained from Ne+U reactions calculated using the hydro- 
dynamic approach with the free streaming freeze-out (up- 
per panel), the Cooper- Frye freeze-out (upper middle panel) 
and the Cooper-Frye freeze-out with the UrQMD afterburner 
(lower middle panel). The lowermost panel shows results from 
default UrQMD calculations without a hydrodynamic phase. 

per nucleon at roughly the speed of sound, which cor- 
responds to an E'kin < 100 MeV, changes the angular 
distribution for all scenarios, in particular for the higher 
beam energies (-Eiab > 1 AGeV) . The high-energetic par- 
ticles predominantly originate from hot and dense regions 
of the head shock which is explicitly excluded by the ap- 
plication of this cut. Therefore, only particles emitted 
from the much cooler and more dilute sideways travel- 
ing Mach shock wave are observed in this analysis. For 
the UrQMD calculations this cut removes most forward 
moving participant nucleons and leaves mostly the fermi 
distributed spectators in the spectrum. 

Plotting dN/d cos 6 instead eliminates the sinusoidal 
shape as can be seen in Fig. [71 Here, the Cooper-Frye 
freeze-out shows a clear peak at very small angles for 
lower beam energies (i?kin < 100 MeV) due to the head 
shock wave and the superimposed thermal distribution 
similar to the one seen in the lowest panel. The results 
from the free streaming freeze-out (uppermost panel) and 
with this the Mach shock wave signal is the less affected 
by this change. 

Figure [S] depicts the dependence of the conical emission 



Figure 6: (Color online) Angular distributions of low energy 
reaction products (iSkin < 100 MeV) of Ne-|-U reactions calcu- 
lated using the hydrodynamic approach with the free stream- 
ing freeze-out (upper panel), the Cooper-Frye freeze-out (up- 
per middle panel) and the Cooper-Frye freeze-out with the 
UrQMD afterburner (lower middle panel). They show emis- 
sion at distinct Mach cone angles in the investigated scenarios 
(cf. Fig. [5]). The maxima of these distributions ^max for the 
different scenarios are depicted in Fig. |8l 



angle ^max on the beam energy, for the different freeze- 
out descriptions considered. The angles are extracted 
from the underlying data of Fig. [6] which includes a low- 
energy cut so that only particles with E'kin < 100 McV 
are considered that originate predominantly from dilute 
and slow fluid cells, i.e. the Mach shock wave. Those 
did not cross the phase transition to the chirally restored 
phase. Therefore, the resulting particle emission at Mach 
angles is probing the purely hadronic phase. It is obvi- 
ous that the Mach wave travels with u ~ Cg through 
the hadronic phase which should result in a rather slight 
change of the emission angle with varying beam energy. 

For the free streaming freeze-out (full triangles), the 
emission angle grows slightly with increasing beam en- 
ergy. For the Cooper-Frye freeze-out disregarding the 
subsequent resonance decays (open squares) however, the 
emission angle is notably shifted towards smaller values 
due to the superimposition of the boosted thermal dis- 
tribution. By employing the UrQMD-afterburner (full 
squares), the computed preferential emission angle gets 
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Figure 7: (Color online) Angular distribution of low energy re- 
action products (Fig. ^ of Ne+U reactions calculated using 
the hydrodynamic approach with the free streaming freeze- 
out (upper panel), the Cooper- Frye freeze-out (upper middle 
panel), the Cooper- Frye freeze-out with the UrQMD after- 
burner (lower middle panel), and default UrQMD (lowermost 
panel), from which the sin(S)-shape of an isotropic thermal 
particle distribution was removed by plotting dN/dcos{6). 



shifted by about A6'iab = +10°. 

Most of the cluster-bound particles from hydrodynam- 
ics with free streaming freeze-out are emitted at Mach 
angles within the range 70° < 6uc < 80° for E'lab > 
1 AGeV [13] . This result corresponds to a cluster flow 
with the speed of sound. 

Comparing the results obtained from the purely hy- 
drodynamic calculations to the results from the hadronic 
transport model calculations without a hydrodynamic 
phase (UrQMD, full circles in Fig. |8]) shows that even 
the pure transport calculation suggests conical emission 
of nucleons, at angles very similar to the expected Mach 
cone angles (note that, for the UrQMD results, the same 
cut in kinetic energy of the nucleons is applied) . While in 
the hydrodynamic picture the conical emission originates 
from a Mach-like wave, binary nucleon-nucleon scatter- 
ing results in the distinct particle emission pattern in 
the transport model which approximates the shape of a 
boosted thermal distribution. 

It is possible to distinguish those scenarios by deter- 
mining the predominant cluster emission angle at non- 
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Figure 8: Beam energy dependence of the emission angle Smax 
for default UrQMD (full circles) and hydrodynamic processes, 
including the free streaming freeze-out (full triangles), the 
Cooper-Frye freeze-out with (full squares) and without (open 
squares) a subsequent UrQMD afterburner. 



vanishing impact parameters. The left panel of Fig. [9] 
shows non-central Ne-|-U collisions at i^iab = 5 AGeV for 
hydrodynamic calculations followed by the free stream- 
ing freeze-out as well as default UrQMD calculations. 
Here, only nucleons emitted within the energy range of 
10 MeV< Skin < 100 MeV are taken into account in 
order to ensure that nucleons originating from the spec- 
tator part of the target nucleus are excluded. In the hy- 
drodynamic scenario, the extracted emission angle stays 
at 0MC — 77° up to an impact parameter of 6 = 4 fm. 
For larger impact parameters, ^lab decreases considerably 
to Slab — 60° for 6 = 9 fm. In this case, the projectile 
hits the target at its periphery (cf. Fig. [2]) and the Mach 
shock wave propagates through the whole of the heavy 
nucleus. UrQMD calculations however, show a different 
behavior. Here, the extracted angle Smax grows with in- 
creasing impact parameter, up to Smax = 90°. This dis- 
tribution pattern corresponds to a fully isotropic particle 
emission. 

Thus, the emission angle for semi-central collisions can 
be used to study the reaction mechanism. While emis- 
sion angles above ^max — 80° favour kinetic collision pro- 
cesses to cause conical emission, emission angles well be- 
low this value indicate collective Mach shock waves in 
nuclear matter. 

The right panel of Fig. [5] depicts the maximum emis- 
sion angle of emitted nucleons in A-l-U collisions with 
£'kin < 100 MeV. Here, the emission angle stays constant 
at 0max — 77° from a up to Ca projectiles. The parti- 
cle number in the low energy bin drops by about 75%, 
and collisions with Ag-projectiles show that there are vir- 
tually no particles with kinetic energies below 100 MeV 
left. The whole dense mater system formed by the two 
colliding nuclei merges completely into a highly accel- 
erated shocked fireball jUj- Compared to these hydro- 
dynamic calculations, the emission angles for A-l-U col- 
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Figure 9: Maximum emission angle Smax of emitted nucleons 
from hydrodynamic (free-streaming freeze-out) and UrQMD 
calculations for non-central Ne+U collisions (left) and for cen- 
tral A-l-U collisions (right) at Biab = 5 AGeV. 



lisions computed with the default UrQMD model are 
larger by roughly 10°. 



IV. SUMMARY 

We studied the emission angles 0iab {d denotes the po- 
lar angle between the beam axis and the flux of matter) 
of reaction products in Ne-l-U collisions in a beam energy 
range from i?iab = 0.5 AGeV to 20 AGeV in a hydrody- 
namic approach with different freeze-out methods as well 
as using the UrQMD model. It is shown that a strong 
head shock wave is created during the hydrodynamic evo- 
lution due to the strong compression of nuclear matter in 
the collision zone. This head shock penetrates the target 
nucleus with a supersonic velocity Vsh and reaches very 
high energy as well as baryon densities which are suffi- 
cient to enter the chirally restored state of matter that is 
described by the EoS. 

Behind the head shock a more dilute Mach cone like 
wave (Mach shock wave) evolves and propagates at slower 
velocity (u ~ Cg) through the hadronic state. 

We show that the energy spectra of emitted nucleons 
peaks at low kinetic energies (see Fig. These slow 
particles (i?kin < 100 MeV) are mostly emitted at Mach 
angles which can be extracted from hydrodynamic calcu- 
lations followed by a free streaming freeze-out (cf. Figs. [5] 
to [7]). Those peaks should be visible even for beam ener- 
gies reached at the FAIR at GSI. These rather small emis- 



sion angles for all investigated beam energies lead to the 
conclusion that Mach-like waves travel in the hadronic 
phase only. Therefore, a crossing of the phase transition 
to the chirally restored phase as it may occur in the head 
shock seems not to be accessible through the emission 
angles of the measured particles. 

Applying however a Cooper-Frye freeze-out (with and 
without the UrQMD afterburner) , these peaks are shifted 
to smaller angles when plotting dN/d9 due to the su- 
perposition with the fully isotropic, but Lorentz-boosted 
thermal particle distribution considered in the Cooper- 
Frye freeze-out. Removing the sinusoidal shape of the 
boosted thermal distribution by considering dN /d cos 0- 
distributions (see again Figs. [S] to [7]), these Mach cone 
peaks can no longer be seen. 

While results from standard UrQMD calculations lead 
to emission angles comparable to those of a hydrody- 
namic calculation followed by a Cooper-Frye freeze-out 
for central collisions (cf. Figs. [5]to[7l), they show a no- 
tably different behavior for non- vanishing impact param- 
eters (see Fig.ini left panel). Thus, the emission angle can 
be used to distinguish the underlying process leading to 
a conical emission pattern of nucleons, checking the ex- 
istence of shock waves in nuclear matter. 

Moreover, we show that for different projectiles from a 
up to Ca a distinct particle emission angle prevails over a 
wide range of impact parameters for both hydrodynamic 
and UrQMD calculations (cf. Fig.[Sl right panel). 
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